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Abstract 

Recent work on stochastic interacting particle systems with two 
particle species (or single-species systems with kinematic constraints) 
has demonstrated the existence of spontaneous symmetry breaking, 
long-range order and phase coexistence in nonequilibrium steady 
states, even if translational invariance is not broken by defects or 
open boundaries. If both particle species are conserved, the temporal 
behaviour is largely unexplored, but first results of current work on 
the transition from the microscopic to the macroscopic scale yield ex- 
act coupled nonlinear hydrodynamic equations and indicate the emer- 
gence of novel types of shock waves which are collective excitations 
stabilized by the flow of microscopic fluctuations. We review the basic 
stationary and dynamic properties of these systems, highlighting the 
role of conservation laws and kinetic constraints for the hydrodynamic 
behaviour, the microscopic origin of domain wall (shock) stability and 
the coarsening dynamics of domains during phase separation. 
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1 Introduction 



1.1 Why we are interested in driven diffusive systems 

The investigation of interacting particle systems far from equilibrium has 
shown that one-dimensional driven diffusive systems with short-range in- 
teractions exhibit a remarkably rich variety of critical phenomena. Unlike 
in thermal equilibrium one observes spontaneous symmetry breaking, long- 
range order and phase coexistence in the steady state if the system evolves 
under certain microscopic kinetic constraints or has more than one conserva- 
tion law. A large body of work has been devoted to microscopic stochastic 
lattice models for driven diffusive systems where classical interacting particles 
move under the action of a random force preferentially in one direction.^ 

The ongoing interest in these systems has many reasons. The most ob- 
vious one is derived from a fundamental task of statistical mechanics, viz. 
the desire to understand the emergence of macroscopic collective properties 
from microscopic interactions, with a view on general features such as inter- 
action range (short- or long-ranged), kinetic constraints or the presence of 
conservation laws. It has turned out that parts of that program can be car- 
ried out to a very satisfactory degree in the simplest case of driven diffusive 
systems of identical conserved particles with hard-core interaction. Despite 
their simplicity, these systems exhibit a rich and rather non-trivial dynamical 
and stationary behaviour. For an exactly solvable paradigmatic model, the 
asymmetric simple exclusion process (ASEP, see below) not only the macro- 
scopic nonlinear hydrodynamics have been derived rigorously PQ |3] but also 
detailed information about universal phenomena, including shock diffusion 
[2|, the microscopic origin of the stability of shocks |1| and the dynamical 
structure function ^ could be obtained in the past decade. It is then natural 
to ask what to expect in the presence of more than one conservation law, i.e., 
in systems with several distinct species of particles). 

Secondly, in the absence of a general framework for studying nonequi- 
librium systems (analogous to the usual principles of equilibrium statistical 
mechanics), one needs to understand coarse-grained dynamical properties not 

^Strictly speaking one should call these models mesoscopic as interactions on scales 
below the particle size are replaced by effective interactions. However, the notion micro- 
scopic has become standard and will be used here. It is justified in relation to a truly 
macroscopic description where particle positions are replaced by coarse-grained density 
fields. 
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only for their own sake, but also in order to predict what stationary states 
these systems evolve into. For one-species systems with open systems this 
has led to a theory of boundary- induced phase transitions which provides a 
general framework for a quantitative description of the steady-state selection 
in driven diffusive systems which are in contact with particle reservoirs at 
their boundary. Unlike in equilibrium, boundary conditions determine the 
bulk behaviour of driven diffusive systems in a decisive fashion which can be 
captured in terms of an extremal principle for the current jHl [Z| . The result- 
ing phase diagram for the nonequilibrium steady state is determined by the 
interplay of localized excitations and shocks. Again it is natural to ask for 
principles of steady-state selection and the resulting phase diagram in sys- 
tems with many species of particles which are characterized by a conserved 
current for each particle species. 

A third motivation for studying these systems stems from numerical ev- 
idence which shows that addressing these questions not only leads to an 
encyclopedic accumulation of knowledge. Rather it was found that there is 
exciting new physics in systems with more than one species of particles, in- 
cluding spontaneous symmetry breaking and phase separation phenomena, 
even in translation invariant systems [313 without the beneficial "assistance" 
of open boundaries or static defects ^U] in facilitating phase transitions. It 
also emerged that similar phenomena may occur in single-species systems 
with one conservation law, provided there are kinetic constraints determined 
by a zero-rate condition somewhat analogous to the zero-temperature con- 
dition for long-range order in equilibrium systems with short-range interac- 
tions. Neither the hydrodynamic behaviour of systems with more than one 
conservation law nor the microscopic conditions for the occurrence of critical 
phenomena are well-understood. 

These are some - and by far not all - fundamental reasons to investigate 
translation invariant driven diffusive systems and to clarify the role of con- 
servation laws and kinetic constraints for their dynamical and steady state 
properties. Other important issues include the calculation of exact stationary 
distributions and large deviation functions, nonequilibrium Yang-Lee-theory 
and phase transitions in systems with absorbing states. Some of these top- 
ics are discussed in a complementary review by Evans ^1]. We also refer 
to (121 Cni UH US] where closely related questions in nonconservative sys- 
tems and higher dimensions are treated and to ^H] and references therein for 
phase transitions in systems with continuous state space. Here we focus on 
the issues of hydrodynamic behaviour, the microscopic origin of domain wall 
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(shock) stability, spontaneous symmetry breaking and coarsening of domains 
in one-dimensional conservative lattice systems, with particular emphasis on 
models with two conservation laws and on exclusion processes with one con- 
servation law in so far they are relevant in the context of the nonequilibrium 
bulk phase transitions that we review here. 

1.2 Basic models 

The asymmetric simple exclusion process jSl E] has become a paradigmatic 
example for a driven diffusive system and has begun to attain a status in 
the study of nonequilibrium systems somewhat similar to the role the Ising 
model plays in equilibrium statistical mechanics. In this stochastic lattice 
gas model each lattice site is occupied by at most one particle. Particles 
hop randomly in continuous time to the right neighboring site with rate D^. 
and to the left with rate Di respectively, provided the target site is empty. 
Otherwise the attempted move is rejected. We present this hopping rule as 
follows: 

AO ^ OA with rate 

OA AO with rate A- (1-1) 

Hopping attempts take place independently with an exponential waiting time 
distribution with mean + (Fig. 1). 

These simple rules specify completely the stochastic bulk dynamics of the 
system.^ For a finite lattice with L sites one has to specify boundary con- 
ditions. Mostly commonly studied are periodic boundary conditions, reflect- 
ing boundaries (hopping confined to a box), and open boundary conditions 
where particles may enter and exit the lattice at the boundary sites 1 and 
L respectively with rates 7,5 (see Fig. 1). By choosing a = piDrXi, 
7 = (1 — pi)DiXi as left boundary rates and (3 = {1 — p2)DrX2, S = 
the open system may be interpreted as being connected to particle reservoirs 
with constant density pi at the left (auxiliary) boundary site and density 
P2 at the right (auxiliary) boundary site L + 1 respectively. The parameters 
Ai^2 are introduced to describe a hopping mechanism between the reservoirs 
and the chain which may differ from the hopping inside the chain. From a 
physics point of view this would correspond to activation energies for enter- 

^For a mathematically precise definition of the process, see I17j. 
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ing (leaving) the system which are different from those for hopping in the 
bulk. 

This is the simplest model that incorporates the basic features of a driven 
diffusive system with short range interactions. The short range interaction 
is taken care of by the hard-core exclusion constraint. The randomness of 
the hopping events models diffusive motion of the free particles outside the 
interaction range. The hopping asymmetry corresponds to a driving force 
that leads to a biased average motion and hence to a macroscopic particle 
current even in the stationary state of the system, except in the case of 
reflecting boundary conditions, where the system reaches an equilibrium state 
[T71 ll8j. Throughout this paper we assume a bias in positive lattice direction. 

Notice that each lattice site can be in two states, either occupied or 
empty, and hence the system can be described in terms of occupation numbers 
Uk = 0,1. Implicit in this description is the absence of any internal degree of 
freedom that particles may possess. Hence all particles are indistinguishable. 
The number of particles is conserved in the bulk, but not at the boundaries 
in the case of the open system. The single bulk conservation law gives rise 
to a current via the lattice continuity equation 

d 

—pk=jk-i-jk (1-2) 

for the expected density Pk = {n-k), averaged over realizations of the stochas- 
tic time evolution and also averaged over different initial distributions. For 
the ASEP the current follows straightforwardly from the definition, 

jk = Dr{nk{l - Uk+i) ) - Di{{l- nk)nk+i ). (1.3) 

For periodic boundary conditions there is a family of stationary distri- 
butions which are Bernoulli product measures with density p, i.e., at each 
given site the probability of finding a particle is given by p, independent of 
the occupation of other sites. According to ()1.3p the stationary current 

j = {Dr-De)p{l-p) (1.4) 

is a nonlinear function of the density with a single maximum at p = 1/2. 

If the ASEP is confined to a box (corresponding to "open boundaries" 
with a = /5 = 7 = 5 = 0) the system evolves into an equilibrium state 
where essentially all N particles form a cluster of macroscopic size ~ N 
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with density p ~ 1 and the current vanishes. The density profile has a non- 
trivial form (deviating significantly from or 1 respectively) only at the left 
edge of the cluster, denoted below as "shock" , "domain wall" or "interface" 
respectively, depending on context. The width of this interface is finite on the 
lattice scale, i.e., microscopic. In the grandcanonical ensemble there are no 
correlations, but the density profile has the form of a hyperbolic tangent ^7] . 
The canonical distribution is more complicated, with correlations within the 
interface region, but has a similar density profile [THj. 

The open ASEP has a intriguing phase diagram with a nonequilibrium 
first order transition at pi = 1 — p2 between a low-density phase with bulk 
density p = pi to a high-density phase with bulk density p = P2- There 
are nonequilibrium second-order transitions from both phases to a maximal 
current phase with p = 1/2, irrespective of the boundary densities in the 
square defined by pi > 1/2, p2 < 1/2. The density profiles are non-trivial 
in all phases 120] • At the first-order transition line one has phase coex- 
istence with domains of densities pi^25 separated by a microscopically sharp 
domain wall (shock), the position of which performs a random walk over the 
whole lattice. The exact solution of the stationary density profiles and a the- 
ory of boundary-induced phase transitions which provides a microscopically 
oriented derivation of the phase diagram is reviewed in detail in The 
bulk densities as a function of boundary densities were obtained by Liggett 
PT] using probabilistic methods. More generally, the theory of boundary- 
induced phase transition which takes into account the flow of fluctuations 
has revealed that the bulk density that simple open driven diffusive systems 
select can be obtained from the extremal principle j2i 

j = max j(p) for pi > p2 

j = min j{p) for pi < p2 (1.5) 

involving the stationary current-density relation that can be obtained e.g. 
from measurements or exact calculations in periodic systems where the den- 
sity is conserved. The boundary densities pi^2 entering (jl.5|) are non-universal 
functions of the rates at which particles enter and leave the system. For 
P2 = the extremal principle was first proposed by Krug jH] on a phe- 
nomenological basis. 

One may relax the exclusion constraint to allow for up to m particles 
on each lattice site. This gives rise to the partial exclusion process j22l IIH] 
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with m + 1 states per site. Also particle systems with nonconserved internal 
degrees of freedom such as velocities in traffic flow models 123 have more 
than one possible state per site, but still obey a single continuity equation 
of the form p.2|) . The form ()1.3p of the current, however, strongly depends 
on the microscopic hopping rules. An important model without exclusion 
is obtained by the following simple mapping of the ASEP: Since the order 
of particles is conserved in the ASEP one may regard particles as sites of 
a new lattice gas system and the number of vacancies Uj between particles 
j, j + 1 as occupation numbers on site j. This gives rise to a special case of 
the zero-range process (ZRP) [211 • The particle hopping rates in the ASEP 
turn into the rates of decreasing the number of particles by one unit in the 
ZRP, with hopping to the right in the ASEP corresponding to hopping to 
the left in ZRP and vice versa (Fig. 1). 

The general homogeneous ZRP allows for hopping of a particle from site 
j with a rate Wn that depends only on the unrestricted occupation num- 
ber Uj. Here we shall consider only nearest neighbour hopping (Fig.2) with 
asymmetry factors p, q to the right and left respectively. Below a critical 
density (which may be inflnite, see Sec. 4), the periodic and inflnite system 
has a family of stationary distributions which are product measures where 
the probability of flnding n particles on a given site is given by [211 l2Zj 

1 1 " 

P, = -, Pn=^^"niM (1-6) 

with the "fugacity" z flxing the mean particle density and the one-site 
nonequilibrium analog 

oo n 

^=E^"niM (1-7) 

n=0 i=l 

of the partition function. The density as a function of z is then given by 

p = z^hiZ. (1.8) 

dz 

According to its deflnition via the lattice continuity equation for the ZRP 
the stationary particle current is given by 

oo 

j = {p- q)Y,PnWn = {p- q)z. (1.9) 

n=l 
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The density dependence of j in the stationary state can be obtained by 
inverting the relation p{z) which is a monotonically increasing function of p 
[Hj. Notice that the radius of convergence depends on the choice of rates Wi 
(see Sec. 4). 

In order to describe a system with two different conserved species A, B 
of identical particles (or alternatively: tagged particles or particles with two 
internal states which do not affect its dynamics) one needs a model where 
each lattice site can be found in at least three different states: empty, or oc- 
cupied by either an A-particle or a 5-particle. The most simple extension of 
the exclusion process that accounts for the possibility of two particle species 
may hence be described by the six hopping rates 
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There is no established name for this generic process and we shall refer to 
it as two-species ASEP. Associated with the two conservation laws there are 
two currents defined by 

JtP^=Jk-i-Jk (1-11) 

^Pf =jf-i-jf- (1-12) 

Notice that in general and depend on both occupation numbers n^, nf 
respectively. Hence one has two coupled lattice continuity equations. The 
stationary distribution of this process and hence the current-density relation 
is known only on certain parameter manifolds, see below. 

The natural order parameter that describes the macroscopic state of the 
system is the particle density of each species. Hence for each conservation law 
there is an associated order parameter. Notice, however, that one may have 
conservation laws that are only indirectly related to particles densities which 
by themselves may not be conserved. E.g. in a reaction-diffusion system 
A+5 —>■ |2H| where A and 5-particles annihilate upon encounter (i.e., react 
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into an inert reaction product), the difference s = — still gives rise to 
a single conservation law with an associated current, even though n"^ and 
are not individually conserved. By interpreting A-particles (i?-particles) as 
carriers of a positive (negative) electrical charge, one could speak in this case 
of charge conservation. In the two-species catalytic reaction A + B ^ B -\- B 
PU] the total density n = + \s conserved. Analogously, one may also 
consider models with two conservation laws, but more than three states per 
site. Examples include e.g. two-lane models jHOl or bricklayer models 

m 

There is no answer to the question to which extent or under which circum- 
stances the existence of internal degrees of freedom matters for the macro- 
scopic properties of driven diffusive systems. Since, however, it is now clear 
that the number of conservation laws is important, we shall categorize the 
models in the following according to this property. 

1.3 Closely related models, not covered in this review 

Equivalence to 2-dimensional equilibrium systems: The time- 
dependent one-dimensional stochastic processes discussed above are equiva- 
lent to two-dimensional equilibrium systems, defined by some vertex model 
|33j . The time evolution is encoded in the transfer matrix of the 2-d model, 
e.g., the discrete-time ASEP with a sublattice parallel update corresponds to 
the six- vertex model [SH 11421 EHl IHZj • For three-states models the construc- 
tion is entirely analogous and leads to higher vertex models. Recently also 
integrable vertex models with an unlimited number of states have been inves- 
tigated [3H]. They correspond to the zero-range representation of the ASEP. 
The stationary distribution of a one-dimensional process with L sites cor- 
responds to the equilibrium state of the associated two-dimensional model, 
defined on a strip of dimension L x oo. Processes defined in continuous 
time are derived from the transfer matrix of the vertex model in the same 
way as one obtains quantum spin-chain Hamiltonians ^ . Hence the Markov 
generator of the stochastic time evolution is equivalent to some (usually non- 
hermitian) one- dimensional quantum Hamiltonian. This opens the tool box 
of condensed-matter physics for the study of stochastic dynamics. The ASEP 
fjl.ip . the two-species ASEP ()1.1U|) on a certain parameter manifold |129j and 
also various single-species jSH HO] and two-species reaction-diffusion models 
EHl 1^ 1121 Hn] correspond to integrable models, for which Bethe ansatz 
and related methods yield exact results on the dynamics of the system. 
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We note that by considering the stationary states of one-dimensional 
non-equihbrium systems as equihbrium distributions of a two-dimensional 
systems the occurrence of long-range order and phase separation becomes 
somewhat less mysterious. However, the important question of how these 
phenomena emerge from the microscopic laws of interaction cannot be 
answered by this formal equivalence. We mention in passing that in another 
mapping the steady-state distribution describes the equilibrium properties 
of a directed polymer in a 2-d random energy landscape 

Higher-dimensional nonequilibrium models: The extension of the 
models discussed above to higher dimensions is technically rather straight- 
forward and obviously of importance. For suitably chosen initial distribution 
some features of shocks may be present in higher-dimensional driven diffusive 
systems, but many of the physical properties of the one-dimensional systems 
discussed below are expected to change dramatically not only because of 
the upper critical dimension = 2 for diffusion (which makes mean-field 
behaviour more likely to describe the systems), but also because of the 
absence of blocking effects due to hard-core exclusion. Moreover, there may 
be phase transitions in the unbiased equilibrium counterparts of the model 
which lead to new phenomena in the biased case. Two-dimensional driven 
diffusive systems are reviewed in detail in [12]. 

Many conservation laws: A natural, but for the scope of this review too 
far-reaching question is the behaviour of particle systems with many con- 
servation laws. Not surprisingly, no general picture has emerged yet. An 
interesting generalization of the ASEP comprises lattice systems with par- 
ticles covering more than one lattice site, but moving only by one site in 
each infinitesimal time step [IHl ESI 113 EHl ESI ■ The stationary distribution 
of these systems has still a simple product structure and remarkably these 
models are also integrable. This includes also polydisperse models with par- 
ticles of different sizes (and hence as many conservation laws). 

Another class of models with many conservation laws arises from assigning 
to each particle its own intrinsic hopping rate. In the zero-range mapping 
one thus obtains a process with site-dependent quenched random hopping 
rates. For rates drawn from some distribution the hydrodynamic behaviour 
has been studied in Refs. jSUlEIlEll- For certain distributions the system 
with asymmetric hopping rates undergoes at some critical density a transition 
to a platoon state, where fast particles are trailing a slow one reminiscent of 
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traffic flow. This transition is a classical analog of Bose- Einstein condensation 
[221 IM| and is quite analogous to the condensation transition in ordered 
systems to be discussed in detail below. A similar model with passing of 
particles has also been studied jSSl EH]- Remarkably, the phase diagram 
of the open system both with j2Z| or without passing jSH] has a structure 
similar to that of the usual one-species exclusion process [T^ 120] which can 
be explained by the theory of boundary- induced phase transitions for systems 
with one conservation law |7l EH] • 



The motivation we gave for studying the systems reviewed here addressed 
general questions of non-equilibrium statistical mechanics, with little 
reference to actual realizations where such processes might play a role. Ap- 
plications are actually numerous and include not only quasi one-dimensional 
settings (e.g. molecular diffusion in nanoporous materials such as zeolites 
[nOl inn E21 inni ISll inni , single- flle diffusion of mesoscoplc colloidal particles 
|66j . or ionic conduction in narrow channels [HZl IHHl ^\ but also - through 
various mappings - two and three-dimensional systems. Of course, the 
basic models (jl.ll) . (jl.lOj) can serve only as very crude approximations for 
any real complex system. However, the universality of critical phenomena 
(dynamical and static), of diffusion, of the emergence of shocks and of 
coarsening allow for the study of fundamental properties of real systems in 
terms of simple toy models. It is not the purpose of this article to review 
such applications in any detail, but some significant results are summarized. 

Tracer diffusion: The simplest way of obtaining a system with two con- 
servation laws consists in considering tagged particles in the usual ASEP. 
Tagged particles (=particles of type B) have the same physical properties as 
usual particles, except that they carry a marker which allows for their iden- 
tification, but does not affect the dynamics. Thus one gets the two-species 
ASEP dEH) with 
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In the unbiased case Dao = -Dqa a single tracer particle in a station- 
ary system of density p is predicted to perform anomalous diffusion with a 
mean square displacement {X^{t) ) oc (l-p)/pv^ [701 [HI 1121 • Recently this 
was confirmed experimentally in the investigation of tracer diffusion in zeo- 
lites jni] using pulsed field gradient NMR and in the study of single-file 
diffusion of colloidal particles [UB] . 

Also a driven tracer particle in an environment of unbiased A-particles 
behaves sub diffusively with a square-root power law for the mean square 
displacement [71]. In the fully driven case ()2.1|) the situation is more 
complex. When averaging over random initial states of the system according 
to the weights given by the stationary distribution, the mean square 
displacement was proved to grow linearly in time with a diffusion coefficient 
D = [Dr — Di){l — p) j^. On the other hand, for fixed initial states 
(averaging only over realizations of the process) the variance is expected to 
grow sub diffusively with power t^^^ [7^1 177j . In a finite system with periodic 
boundaries the variance in the number of hops made in the totally asym- 
metric process {Dqa = 0) has been calculated exactly in the infinite-time 
limit [7H1 Eni and been found to decrease asymptotically oc 1/ \/L in system 
size. This is to be expected from dynamical scaling with the well-known 
dynamical exponent z = 3/2 of the asymmetric exclusion process jSUl IHT]. 

Shock tracking: The ASEP exhibits shocks which on a macroscopic level 
appear as stable moving discontinuities in the density profile. It is of great 
interest to understand the microscopic structure of shocks, i.e., the density 
profile and correlations on the microscopic lattice scale which for a real system 
is the analogue of intermolecular distances. The fundamental question is 
whether the density changes quickly over molecular length scales or much 
slower (but still abruptly on macroscopic scales). Trying to answer this 
question leads to the problem of defining a microscopic shock position in a 
given realization of the process. This can be accomplished by introducing 
a second class particle, i.e., a particle that moves w.r.t. vacancies like an 
ordinary (first class) particle, but behaves like a vacancy w.r.t. to ordinary 
particles [H2]- This leads to rates (jl.lUp with 
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By studying the motion of a single second-class particle one finds the mean 
shock velocity 

Vs = (2.3) 

Pi- P2 

for a shock jumping from a density pi to p2 with stationary current ji^2 
in each domain. The expression ()2.3|l may be deduced directly from mass 
conservation. For fixed initial states the variance of the shock position is 
subdiffusive with power law t^^^ [76 , while with averaging over random initial 
states with stationary weights at different densities pi^2 to the right and left of 
the shock (i.e., starting the system from a shock measure) one finds ordinary 
diffusion. The diffusion coefficient was conjectured |83^ and subsequently 
proved jHlj to be given by 

D. = iliiii. (2.4) 

2 Pi - P2 

More detailed information about the microscopic structure of the shock has 
been proved by a variety of methods, for a review see and for more recent 
work [SHllHnilHZl and references therein. It has been established that a shock is 
truely microscopic in the sense that a rapid increase of the density of particles 
occurs on the lattice scale. Loosely speaking one may say that the shock 
performs a random walk with drift velocity (j2.3|) and diffusion coefficient 
(j2.4j) . Some details concerning the structure of a shock are reviewed below. 

In a region of smooth variation of the density the second class particle 
allows for tracking localized perturbations of the density j3 EH] • The mean 
velocity of the second class particle is given by the collective velocity 

^c=^j(p) (2.5) 

of a density perturbation. Its mean-square displacement grows superdiffu- 
sively with power law 

{X\t)) - {X{t)f oct^/\ (2.6) 

where the divergent effective diffusion coefficient A crosses over in a finite 
system to jHO] 

A ~ (2.7) 

which is expected from dynamical scaling with dynamical exponent z = 3/2. 
As a toy model for econophysics the position of the second class particle 
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marks the price of an asset on the price axis in a hmit order market, the 
first-class particles represent bid and ask prices jSHl- The rigorous result 
()2.6|) corresponds to a Hurst exponent H = 2/3 for fiuctuations which 
compares well with the empirically observed value H ^ 0.6 for intermediate 
time ranges. Introducing annihilation and creation of particles (cancellation 
and renewal of orders) leads to the Gaussian value H = 1/2 after some 
crossover time, also in agreement with empirical findings. 

Traffic flow: The occurrence and microscopic nature of shocks in the ASEP 
is reminiscent of traffic jams in vehicular traffic. Indeed, traffic flow may be 
regarded as a driven diffusive system [211 123; albeit with a non-conserved 
internal degree of freedom, viz. the speed of cars which is dynamically de- 
termined by the competition of the desire to move at an optimal high speed 
and the necessity to keep a velocity-dependent minimal safety headway to 
the next car. At low densities the mean distance between cars is larger than 
the required safety headway and essentially all cars move at their optimal 
speed, with some fluctuations. At high densities cars have a mean distance 
below the safety headway corresponding to the optimal speed which leads to 
a slowing down of the traffic. As a result the mean current as a function of 
density p (known as fundamental diagram in the traffic literature) has a max- 
imum like the exclusion process, albeit with a much broader distribution of 
the current. In the first measurement of traffic flow in 1935 by Greenshields 
[DTj the measured mean flow of cars was approximated by the expression 
()1.4|) j oc p{l — p) that the ASEP yields. More sophisticated models that 
provide a much better description of real traffic data have been developed 
in the past decade, starting with the Nagel-Schreckenberg model [121 which 
contains the ASEP as a simple limiting case. Yet some fundamental features 
of the ASEP survive in the more complicated Nagel-Schreckenberg model. 
The theory of boundary-induced phase transitions developed in fJl ^] which 
predicts the stationary phase diagram of a single-species system with open 
boundaries in terms of the extremal principle ()1.5|) explains quantitatively 
the phase diagram of the Nagel-Schreckenberg model in terms of effective 
boundary densities and is also consistent with measurements of real highway 
traffic j93j. 

A description of traffic flow with a lattice model with one conservation 
law corresponds to modelling cars which all have the same intrinsic optimal 
speed - a rather crude approximation if describing mixed traffic of cars and 
trucks is envisaged. This naturally leads as a next approximation step to 
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a two-species description with "fast" and "slow" particles respectively. A 
model with two conservation laws arises also from the study of two-way 
traffic flow with interaction between lanes, but no exchange of particles in 
each lane. Taking the exclusion process as the simplest possible model one 
arrives at the model (jl.lOp with suitably chosen rates [Hll ESI EH] , see below. 
A mixture of cars with individual intrinsic speeds leads to the disordered 
hopping models mentioned above. 

Biophysics: The ASEP with open boundaries was first developed as a sim- 
ple model for describing the kinetics of protein synthesis [HZl HH] • Here parti- 
cles are ribosomes moving along the codons of a messenger RNA. Each codon 
corresponds to a specific aminoacid which the ribosome uses to assemble a 
protein. When such a step is completed the ribosome moves on to the next 
codon and continues with the addition of next aminoacid to the growing pro- 
tein molecule. The injection of particles at one end marks the initialization of 
the process, the absorption at the other boundary describes the release of the 
ribosome. The shock known from the exclusion process corresponds to a traf- 
fic jam of ribosomes which explains an experimentally observed slowing down 
of the ribosomes as they approach the terminal point of the m-RNA where 
they are released after completion of the protein synthesis |1H1 UHl 1^^! HUUj • 

Very recent work has shown that in another biological setting exclusion 
particles may describe molecular motors such as kinesins moving along micro- 
tubuli or actin filaments in a cell |1(JHI1(J2] . Due attachment and detachment 
during the motion a description with non-conservative dynamics where parti- 
cles are annihilated and created also in the bulk with a small rate is required. 
This leads to the model of Ref. jBHl with open boundaries which yields in- 
teresting new phenomena |lU3j . Oppositely moving molecular motors give 
rise to a two-species exclusion model, in analogy to two-way traffic fiow with 
interaction between lanes. 

A two-species exclusion process has also been introduced to describe the 
motion of ants along ant trails |lU4j . While crawling along a trail, the ants 
- modelled by A-particles hopping along a lattice - produce pheromones 
(i?-particles) which serve as a marker of the traversed path for other ants 
which again produce pheromones for subsequent ants. This is necessary 
to stabilize the trail as the pheromones evaporate after some time. The 
pheromones are modelled as an immobile particle species which is deposited 
when a hopping event has taken place and which disappears with some 
evaporation rate. The analogy of the fiow of ants to traffic fiow has been 
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pointed out in Ref. |l(J5j who measured the flow rate versus the density of 
ants, i.e., the current density relation. The numerical results obtained from 
the two-species ant trail model yield qualitatively similar results |lU4j . Es- 
sentially the same model (with different parameter values and update rules) 
has been introduced as a "bus route" model where one observes bunching of 
particles (= "buses") as they travel along lattice ("bus stops") and pick up 
passengers |lU6j . Bunching of real buses appears to occur on services which 
do not run according to fixed schedules, but which stop according to demand. 

Polymer dynamics: In polymer networks such as rubber gum or gels, in 
polymer melts or in dense solutions of macromolecules such as DNA different 
polymer strands form a complicated topological structure of entanglements 
somewhat reminiscent of a large portion of spaghetti. The entanglements 
severely restrict the dynamical degrees of freedom of the polymer chains. In 
the framework of the celebrated reptation theory |107j developed by Doi, 
Edwards, and de Gennes |1U81 llU9j the motion of an individual polymer is 
viewed as being confined by a hypothetical tube which models the collec- 
tive effect of all entanglements of the neighbouring polymer chains. In an 
uncrosshnked melt or solution the tube is open at both ends, since at the 
end points the motion of polymer segments transverse to its own contour is 
not restricted by topological constraints. This picture results in a snakelike 
one-dimensional effective dynamics of polymer segments along the tube, with 
extra orientational degrees of freedom only at its ends. 

In a lattice model of Rubinstein [IlUj the reptation dynamics is modeled 
by the symmetric exclusion process with open boundaries which de- 

scribe the extra end point degrees of freedom. With this model exact results 
for the relaxation of the contour and contour length fluctuations have been 
obtained. Recent experiments on the dynamics of single entangled DNA- 
molecules in dense solution confirm the findings |11H I112j . 

Duke |113j extended the model to allow for tracking the spatial orientation 
of the tube rather than only its length. This was done in order to introduce 
a reference axis for describing gel electrophoresis, i.e., the separation of poly- 
mer fragments according to their length L. By applying an electric field of 
strength E (the direction of which is the reference axis) a charged polymer 
is expected to move through a gel matrix (which provides an entanglement 
network) according to the rules of reptation. However, standard reptation 
theory does not allow for a prediction of the drift velocity v beyond the linear 
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response regime of small fields or very long polymers where 

V oc DEL. (2.8) 

Here D is the diffusion coefficient of an entangled polymer, predicted by 
reptation theory to scale 

D oc l/L^ (2.9) 

with length. The extended Rubinstein-Duke model is an asymmetric three- 
states exclusion process (jl.lOp with Dab = Dba = and open boundaries. 
Exact and rigorous results |114[ I115j confirm the predictions ()2.8p . ()2.9|) . 
Moreover, simulations at high fields yield the drift velocity in the non-linear 
regime which are in good agreement with experimental data |117j . At 
sufficiently high fields the model exhibits spontaneous symmetry breaking in 
the orientation of the polymer chain |118j . The asymptotic behaviour ()2.9p of 
the diffusion coefficient has also been proved to remain valid in the presence 
of quenched kinematic disorder which is described a generalized Rubinstein 
model with many conservation laws where particles have their individual 
hopping rates |119j . 

A long-standing mystery in reptation theory has been the asymptotic 
behaviour of the viscosity 77 of a polymer melt which is expected to scale 
asymptotically [TMITnn] 

r](xL^. (2.10) 

However, experiments consistently give higher value ~ 3.4 of the scaling 
exponent. Doi had suggested this to be a finite-size effect due to tube- length 
fluctuations |lU8j . That tube- length fluctuations lead to an increased 
effective exponent could be confirmed by a careful numerical analysis of the 
Rubinstein-Duke model |12(Jj . Also details of the end-segment dynamics 
were shown to have significant non-universal impact on finite-size behaviour 
of the viscosity and the diffusion coefficient |121j . 

Spin relaxation: By interpreting occupation numbers as classical spin vari- 
ables Sfc = 1 — 2nk the ASEP describes biased Kawasaki spin-exchange dy- 
namics |122j for the one-dimensional Ising model at infinite temperature. At 
finite temperature, biased Kawasaki dynamics correspond to an exclusion 
process 

XAOY XOAY with rate D^^ 

XOAY ^ XAOY with rate Df^. (2.11) 
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with next-nearest neighbour interaction depending on the occupation X, Y. 
With some constraints on the rates jHZj this model has a stationary distri- 
bution 

P^n) = le-'3(^(^)+^^) (2.12) 

Zj 

which is an Ising measure with energy E = —J ^i^rikUk+i and a magnetic 
field /i which plays the role of a chemical potential in the lattice gas interpre- 
tation of the variables. At sufficiently low temperatures the particle (=spin) 
current becomes a non-convex function of the density (=magnetization) 
and unusual phenomena such as splitting of shock fronts which separate 
regions of different density can be observed [3 I123j . In spin language the 
three-states model ()1.1U|) describes dynamics of a classical spin-1 system. 

Interface growth: It was already realized in the 1980ies that the ASEP 
describes the dynamics of a fluctuating interface by considering the spin 
variables as local discrete slopes of an interface on a two-dimensional sub- 
strate |124l 1125] (Fig. 2). Hopping of a particle to the right between sites 
k,k + 1 corresponds to the random deposition of a particle on site k the 
dual growth lattice, hopping to the left to an evaporation (Fig. 2). One 
thus obtains a growth model in the universality class of the one- dimensional 
KPZ equation |126j . reviewed in W^^. We stress that the mapping is not 
one-to-one. Since in the exclusion presentation only the local slopes enter, 
the information about the actual height of the interface gets lost. One can 
keep track of the height by introducing an extra random variable for the 
local height at some reference point /cq, which is increased (decreased) by 2 
units whenever a particle hops across the bond /cq, k^ + l to the right (left). 
The steady state current of the ASEP then gives the average growth velocity, 
while fluctuations of the current measure fluctuations in the local interface 
height. The extension of this mapping to the generalized exclusion process 
()1.10|) is obvious, one obtains a system where local height differences may 
take values 0, ±1. Some growth dynamics considered below have the prop- 
erty that particles cannot be chipped off a complete layer. This corresponds 
to a hidden conservation law, which cannot be expressed in terms of particle 
occupation numbers alone. 
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3 Steady states and hydrodynamic limit 



3.1 Steady states for driven diffusive systems 

As has become clear above, the stationary behaviour, i.e., the state the sys- 
tem evolves into, is the first question to be addressed in the investigation of 
driven diffusive systems.^ In this article we are concerned with the behaviour 
of translational invariant systems, defined either on a finite lattice with peri- 
odic boundary conditions or on the infinite integer lattice 2. Consequently 
we shall investigate stationary distributions which are either translation in- 
variant or where translation invariance is spontaneously broken. 

It is important to bear in mind that entirely different dynamics may have 
the same stationary distribution. Indeed, for any given distribution one may 
always construct some equilibrium dynamics (obeying reversibility) using the 
principle of detailed balance. Moreover, a strongly nonequilibrium system 
may have the stationary distribution of some equilibrium model, an example 
being the KLS model 1)2.111) . Therefore, equality of stationary ensembles for 
different systems has no implications whatsoever on the dynamical proper- 
ties of these models. We also remark that simple dynamical rules may result 
in stationary distributions with a complicated structure and long-range cor- 
relations (see below), while complicated dynamical rules may very well lead 
to simple stationary distributions. Some stationary distributions for one- 
species models have been reviewed in the introduction, here we focus on the 
two-species ASEP. 

Fortunately, not only the single-species ASEP but also many two-species 
stochastic particle systems of interest have simple stationary distributions, 
the simplest being product measures with stationary probabilities of the form 

p^gM^ivW^^ (3.1) 

Here A^"^'^ = J2k ^t'^ conserved total particle numbers of each species 

(for one-component systems one has nf = 0) and /i"^'^ are the corresponding 

■'This is formally analogous to investigating the equilibrium behaviour of a many-body 
system and hence in the mathematical literature stationary states are often referred to 
as equilibrium states, even though the presence of macroscopic currents prevents the ap- 
plicability of the usual notions of equilibrium statistical mechanics such as reversibility 
and detailed balance. We note, however, that as long as only a stationary distribution 
is concerned - without reference to the stochastic dynamics for which the distribution is 
stationary - it is convenient to use notions borrowed from equilibrium statistical physics 
such as partition function or canonical/grandcanonical ensembles respectively. 
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chemical potentials. For fixed A^"^, A^'^ all configurations are equally likely. 
In equilibrium, such measures correspond to non-interacting systems. 

By writing the master equation for the stochastic dynamics in the quan- 
tum Hamiltonian formalism ||4j it is straightforward to determine what dy- 
namics have stationary product distributions. In this formalism a product 
measure is represented by a tensor product vector \P) = \p)^^. Each factor 
I p ) in the tensor product has as components the probabilities of finding a 
lattice in a given state. Hopping events between sites k,k + 1 are generated 
by a local stochastic matrix hk acting nontrivially only on the terms k,k + 1 
in the tensor product. The full time evolution is generated by the stochastic 
Hamiltonian H = J2k and the stationarity condition for | P ) reads 

H\P) = 0. (3.2) 

Because of translational invariance a stationary product measure therefore 
satisfies the relation 

hk\P) = {dk+i-dk)\P) (3.3) 

with an arbitrary matrix dk acting nontrivially only on site k. This relation 
is usually very easy to verify. A similar approach can be chosen for more 
complicated measures, e.g., Ising measures with stationary probabilities of 
the form 

p oc e-^Efe(^^^"^ti+-^''''«i+^''''«+i)+A'^^-'+'^''^''. (3.4) 

Given these measures the current can be calculated exactly as a function 
of the densities p^, via the invertible relationship between the chemical 
potentials and the densities. For the 2-states ASEP ()1.10|) one has a product 
measure on the parameter manifold defined by 

DoA - Dao + Dbo - Dob + Dab -Dba = (3.5) 

For a i^'-species ASEP with rates Dxy for the hopping process XY YX 
there are {K — 1) {K — 2) /2 conditions for the existence of a product measure 
[121 

Dox - Dxo + Dyo - DoY + Dxy - Dyx = 0. (3.6) 

In the quantum Hamiltonian formalism the stationary distribution is the 
ground state vector of the associated quantum spin chain. For integrable 
models jl28t I129j one may use the Bethe-ansatz and symmetry properties 
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[TH] for the explicit construction of stationary states which are not simple 
product measures. Using the Bethe ansatz has not been attempted yet for 
this class of models. A more popular method is the application of the matrix 
product ansatz, reviewed in |13Uj . In this approach one defines a product 
measure with matrix entries Dm rather than c-numbers as stationary weights 
for finding a given site in state m. The matrices Dm to together with a set 
of auxiliary matrices have to satisfy algebraic relations which are 

obtained from requiring the matrix product state to satisfy the stationarity 
condition ()3.2|) . This leads to algebras with quadratic relations |2D1 EHl II2Z1 
una .135, 136, 137, 138J. 

The matrix product construction is equivalent to writing the stationary 
distribution of the lattice gas in terms of a transfer matrix C = Y.m=o of 
some n-states equilibrium system, determined by the representation of the 
algebra, in particular its dimension n which may be finite or infinite. It is 
clear that finite-dimensional representations correspond to stationary states 
with exponentially decaying correlations unless the largest eigenvalue of C 
is degenerate. Usual product measures (complete absence of correlations) 
correspond to one- dimensional representations. For a detailed review, see |1H 
I130j . The approach can also be extended to describe the full time evolution 
and hence yield time-dependent probabilities of the system evolving from 
some nonstationary initial distribution |13H 113!^ I133j . Popkov et al. have 
identified the parameter manifold for which the dynamics of the 2-species 
ASEP can be solved using the dynamical matrix product ansatz |129j . Such 
models are all integrable in the sense of being associated with an integrable 
vertex model. The more general models of Ref. |127t I138j for which only the 
stationary distribution can be constructed with matrix products include also 
non-integrable models. 

3.2 Steady states with one 5-particle 

A series of intriguing results have been obtained for two-species systems 
(I1.1U|) with just one particle of type B. Conditioning on having a second 
class particle at some given site and calculating the probability of finding a 
first-class particle at distance r yields the density profile as seen from the 
shock position, defined by the position of the second-class particle. The den- 
sity approaches its asymptotic shock densities pi,2 at an exponential rate, 
given in a non-trivial way by the hopping asymmetry D^/Di and the densi- 
ties pi^2 I139j . For a special value of the asymmetry one has Bernoulli 
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measures with densities pi^2 to the left and right respectively. This is the 
result of a g-deformed S't/(2)-symmetry of the Heisenberg quantum Hamil- 
tonian that generates the time evolution of the process. For this value of 
the asymmetry (or, equivalently, arbitrary asymmetry, but special density 
P2) the time evolution of the shock measure has been calculated exactly both 
for the continuous-time ASEP and for a discrete-time variant jnZl. The 
shock position performs a lattice random walk with rates given by the cur- 
rents and densities in the two branches of the shock. Two consecutive shocks 
with densities pi,2,3 which can be defined by two second-class particles form 
a usual bound state with finite mean distance and exponentially decaying 
distance distribution if a condition on /92,3 originating in the g-deformed 
SU (2)-symmetry is met |14Uj . Generically two second-class particles form 
a weak bound state with infinite mean distance and algebraically decaying 
probability p{r) oc r~^/^ of being a distance r apart [H^ I139j . 

Using the algebra arising from the stationary matrix product ansatz 
Mallick has studied the two-species ASEP p.lOj) with rates 



All other rates are zero. For a = 1, (3 = the "impurity" particle B cor- 
responds to a tracer particle, for a = /? = 1 it is a second-class particle. 
This model describes "cars" (A-particles) and "trucks" (-B-particles), with a 
passing rate f3. For a < 1 and (3 < a a single truck acts like an impurity, 
hindering the motion of cars. The current of A-particles, the velocity of the 
impurity and the density profile as seen from the impurity have been calcu- 
lated exactly [HI]. The system with a single impurity exhibits an interesting 
phase diagram as a function of the hopping rates. In one of the phases the 
system develops a stationary shock for sufficiently large density p > pc, anal- 
ogous to a traffic jam building up behind a slow vehicle. One has coexistence 
of a low-density domain and a high-density domain, separated by a domain 
wall (Fig. 3). The diffusion constant of the impurity has also been calculated 
exactly pmiTlTj. 

Lee et al. [0^] considered the model (jl.iup with rates 



corresponding to oppositely moving particles (slow-moving "trucks" for 7 < 
1) which interact upon encounter. In the presence of a single truck (impurity) 
the average speed of cars (and hence the current), the speed of the truck and 



Dao = 1, Dbo = a, Dab = 



(3.7) 



Dao = 1, Dob = 7, Dab = l/P 



(3.8) 
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the density profile of cars lias been calculated using again the same stationary 
three-species algebra [HSl- One obtains a phase diagram with a transition to 
a jammed phase at a critical density Pc = (Fig- 3), with the remarkable 
property of having the same statistical properties as a deterministic ASEP 
with a fixed impurity |142j . In a finite system of site L the position of the 
(microscopically sharp) domain wall fluctuates over a region of length oc y/L. 
In a system with two trucks they form a weak bound state in the traffic 
jam phase. Notice that the model ()3.8p is equivalent to ()3.7|) by exchanging 
S 0. However, in this mapping a single truck corresponds to a single 
vacancy, a scenario not studied in W^ . 

Arndt et al. |143j introduced a model of type (j3.8p with 

Dao = A, Dob = A, Dab = q, Dba = 1 (3.9) 

For finite densities of both particle species this model is reviewed below. 
Jafarpour [2^] considered the presence of a single S-particle and calculated 
exactly for A = 1 the speed of A-particles and the impurity and the density 
profile of A particles. As in the model of Lee et al. there is a phase transition 
from a free-flowing to a jammed phase, here at a critical density pc = q/2 
(Fig. 3). 

We remark that the jamming transition seen in these three versions of 
the two-species ASEP may be regarded as a kind of condensation transition 
where a finite fraction of A-particles condenses into a macroscopic block 
trailing the moving impurity. A similar transition occurs also in the usual 
ASEP with a fixed blockage |142H144]ll45j . corresponding to an immobile B- 
particle and passing event ABO OB A rather than AB BA. Hence single 
particles in a system with two conservation laws play a role somewhat similar 
to local inhomogeneities in a system with one conservation law. Notice that 
for a fixed blockage and for the model of Mallick the jammed phase exists 
between two critical densities p~,p^. Lee et al. and Jafarpour resp. report 
only a lower critical density, below which the system is in the free-flow phase. 

3.3 Hydrodynamic limit for finite densities 

The link between the stationary state and the dynamics is established by the 
continuity equation which relates the change in the local order parameters 
to the currents and other equations for the dynamics of the various non- 
conserved internal degrees of freedom. In order to obtain information about 
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the dynamics on a course-grained scale one assumes local stationarity and 
a sufficiently smooth behaviour of the local order parameter. For driven 
diffusive systems one may then investigate the dynamics on the Euler scale, 
i.e., in the scaling limit where the lattice spacing a and time scale r are sent 
to such that a/r remains constant. 

By expressing the correlation functions that enter the currents completely 
in terms of the density one thus obtains a partial differential equation 

d d , ^ , ^ 

for systems with one conservation law, and system of equations 

- -h'^P-'P") (3-11) 

for models with two conservation laws. Using the theory of partial differential 
equations this yields coarse-grained information about the time-evolution 
of the lattice model. The average occupation of the local density is thus 
described in terms of a deterministic evolution of a coarse-grained density 
profile. 

From these introductory remarks it has become clear that knowing the 
stationary currents exactly as functions of the density is crucial for calculating 
the density profile. A mean-field approximation which neglects correlations 
between different lattice sites is bound to give wrong quantitative results 
unless the stationary distribution happens to be characterized by the ab- 
sence of correlations. For systems with a single conservation law mean field 
approximations may yield qualitatively correct behaviour if the mean-field 
current reproduces local extrema and infiection points of the true current. 
However, we wish to stress that even short-ranged correlations resulting from 
Ising-type stationary distributions ()2.12|) may lead to qualitatively wrong be- 
haviour of the current- density relation. An example is the KLS-model (|2.1H) 
where the exact current has a local minimum between two symmetric max- 
ima [1.461 whereas the mean-field approximation yields the current-density 
relation ()1.4|) of the ASEP with a single maximum. Hence neglecting short- 
ranged correlations may yield not only quantitatively but even qualitatively 
wrong predictions for the most basic dynamical properties of the system, viz. 
the stationary current and the coarse-grained time evolution of the density 
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profile. To conclude, one needs to know which generic features of a current 
density relation determine the qualitative behaviour of the solution of the 
hydrodynamic equation and it must be verified that a mean-field approxi- 
mation reproduces those features. Otherwise a mean field treatment of the 
continuity equation yields no information about the dynamics of the sys- 
tem. In some cases approximating an unknown measure by an Ising measure 
(j2.12|) . (j3.4jl with short-ranged correlations rather than by a simple product 
measure (simple mean field) without correlations may bring improvement. 
In the literature this improved kind of approximation scheme is sometimes 
called cluster approximation. In the following we assume that at least the 
stationary current (if not the full measure) is known exactly. 

Nonhnear equations of the form (|3.1(J|) , (|3.11|) are known to possess singu- 
larities which do not allow for a unique solution of the initial value problem. 
Almost all initial configurations will develop shock discontinuities where the 
density jumps from one value to another even if the initial state was smooth. 
This raises the question of the microscopic properties of the macroscopic so- 
lution. Moreover, ambiguities exist even if the original particle problem has 
a unique stationary state into which the system evolves for all initial states. 
Thus one is faced with the second problem of selecting the physical solution 
of the hydrodynamic equation (Fig. 4). 

For a class of systems with one conservation law the transition from the 
stochastic lattice dynamics to the hydrodynamic equation (|3.1(J|) is math- 
ematically well-understood f£i7\ and includes also the treatment of shock 
discontinuities, for a broader overview see p. Here we give a more physics 
oriented account inspired by the desire to derive macroscopic phenomena 
such as shocks from microscopic behaviour, viz. the flow of localized (micro- 
scopically) perturbations inside a stationary region. It has turned out that 
such an approach, originally developed for systems with one conserved den- 
sity and no internal degrees of freedom, also works for two-species problems, 
if suitably generalized. 

To keep the discussion simple we restrict the review to a system with 
a convex current-density relation such as found in the ASEP (|1.4p . In the 
hydrodynamic limit (|3.1Up one obtains for the ASEP the well-known inviscid 
Burgers equation |148j 

= -g^iDr - A)p(l -p) = -{Dr - A)(l - 2p)—p. (3.13) 
It is well-known that an upward shock in the ASEP constitutes a stable 
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shock whereas the downward shock dissolves into a rarefaction wave. This 
result can be obtained by using the method of characteristics which trace the 
motion in space-time of points of constant density. One introduces a scaling 
variable u = x/t to obtain from ()3.13|) 

M = (D, - A)(l-2p). (3.14) 

This is a weak solution for an initial profile with a down-shock at a; = 
(Fig. 4). 

This solution may be obtained in a different way by starting from the 
lattice continuity equation (jl.2|) with current 

jk = {Dr - Di)pk{l - Pk+i)- (3.15) 

obtained from (jl.3j) by neglecting correlations at all times. The picture un- 
derlying this approximation is the assumption of local stationary where all 
correlations are sufficiently small. With a Taylor expansion in the lattice 
constant a — one arrives at the Burgers equation 

|p=-|(Z,,-AWl-rt + .|jP (3.16) 

with an infinitesimal viscosity v o: a. The Burgers equation is integrable 
using the Hopf-Cole transformation p = K,dx\TLW. This leads to a standard 
linear diffusion equation for w which has a unique solution for any initial 
profile. Taking the limit a — in the solution one recovers the physical 
solutions described above (shock and rarefaction wave) which are realized by 
the ASEP. 

In a different, but for one-species models equivalent approach for selecting 
the physical solution one defines an entropy function associated with the 
conservation law ()3.10|) . The entropy solution yields the physical solution 
corresponding to the ASEP, for a review, see . With this approach one may 
also consider the ASEP with a localized defect |149j which has a different 
solution for a downward shock and which also produces a shock inside a 
domain with constant density [ISUj . 

In order to obtain a physical microscopic picture of how these solutions 
emerge on a macroscopic scale we study the dynamics of localized pertur- 
bations in a homogeneous stationary environment j3]. The time evolution 
of such a perturbation on the lattice scale can be probed by examining the 
dynamical structure function 

S{k,t) = {nk{t)no{0))-p' (3.17) 



26 



which measures the density relaxation of a local perturbation in the station- 
ary state of uniform density p. Generally, the width of such a perturbation 
at t = is of the order of the bulk correlation length. The center-of-mass 
velocity of the perturbation is given by the collective velocity 

^c=^j(p) (3.18) 

One may derive this relation from the shock velocity 

Ji - 32 



Pi- P2 



(3.19) 



by taking the limit pi — > p2 of the asymptotic densities of the shock. Notice 
that Vc changes sign at local extrema of the current density relation. We 
stress that the expression (j3.19|) follows from mass conservation and hence no 
specific assumptions on the nature of the microscopic dynamics are involved. 
Hence also ()3.18|) is expected to be generally valid. The only assumption are 
sufficiently rapidly decaying bulk correlations in the steady state as otherwise 
the microscopic definition of a shock position becomes questionable, since 
the shock position cannot be defined on a scale below the bulk correlation 
length. A direct derivation of ()3.18|) from the dynamical structure function 
which uses only translational invariance, the conservation law and decay of 
correlations is outlined in jl]. 

The shock velocity for the ASEP follows from (jl.4j) . (j3.19|) and one finds 

Vs = iDr-Di){l-p^-p2). (3.20) 

This yields the collective velocity 

v, = {Dr.-Di)il-2p). (3.21) 

The origin of the physical solution of the macroscopic time evolution can 
now be explained from a microscopic viewpoint by imagining that by a small 
fiuctuation a certain amount of mass detaches from the shock and forms a 
perturbation at a small distance from the shock position. On average this 
fiuctuation will then travel with speed (j3.21|) where p is to be taken as either 
pi or p2, depending on whether the fiuctuation had originally moved to the 
left or right of the shock. Eq. ()3.21|) shows in the case of an upward shock 
that for all shock densities pi,2 

vi^'^ > Vs > vf\ (3.22) 
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Hence in the moving reference frame of the shock the excess mass drifts back 
to the position of the shock and hence stabihzes it. 

On the other hand, in the case of an initial downward shock the fluctuat- 
ing excess mass moves on average away from the shock. Therefore this shock 
is not stable against fluctuations, in the course of time the shock smears out 
and develops into a rarefaction wave. In order to predict the macroscopic 
shape of the rarefaction wave, we assume an initial configuration with a shock 
with densities pi and p2 which is composed of many infinitesimal subsequent 
shocks at various levels of intermediate densities (Fig. 3). Neither of these 
shocks is stable, but each slowly dissolving shock at density p moves with a 
speed Vc- Hence we conclude that on the Euler scale (where the spread of a 
perturbation (|2.6p and hence of the increasing width of the unstable shock 
is scaled to zero) points of constant density p generally move with speed 
Vc- From this observation the explicit form of the rarefaction wave can be 
constructed. 

Clearly, this is not a rigorous argument. Support for this picture comes 
from the hydrodynamic limit. The collective velocity is then nothing but the 
speed of the characteristics of the corresponding hydrodynamical equation 
dtp = —dxj resulting from the continuum limit of the lattice continuity equa- 
tion ()1.2|) . In this limit, the criterion ()3.22j) becomes the defining property 
of a shock discontinuity |151j . It asserts that the characteristics are moving 
into the shock. Otherwise, the characteristics yield the rarefaction wave, as 
rationalized above. For current- density relations which are not globally con- 
vex one decomposes a single shock into subsequent small shocks and then 
applies ()3.22|) to these minishocks in order to decide on stability. By taking 
the limit of infinitesimal shocks one recovers in this way the scaling solution 
of the hydrodynamical equation obtained from the method of characteristics 

For systems ()3.1H) with two or more conservation laws there is no well- 
established mathematical theory for the selection of the physical solution 
in the corresponding lattice gas. In recent work 32j Toth and Valko have 
obtained rigorous results by making use of Yau's relative entropy method 
|153j which essentially proves that a product measure with time-dependent 
local densities evolving according to the solution of the hydrodynamic equa- 
tion converges to the true measure in the sense of relative entropy of the 
two measures. This approach works for systems with a stationary product 
measure until a shock has formed, provided some particular identities hold 
which relate the macroscopic fluxes in the hydrodynamic pde. These iden- 
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titles are reminiscent of the Onsager's reciprocity relations. The systems 
studied in [32] are models with generically more states than conservation 
laws. They include a family of 2-species ASEP's with parameters satisfying 
the relation ()3.5p guaranteeing the existence of a stationary product measure 
and which have a natural interpretation as growth models. A two-species 
zero-range process has been studied by Grosskinsky and the corresponding 
hydrodynamic equations have been established |154j . 

The rigorous approach of |22] is rather powerful but in its current state 
fails as soon as shocks develop. On the other hand, shock waves and special 
rarefaction waves have been analysed nonrigorously from a microscopic view- 
point by studying the flow of perturbations and correspondingly extending 
the physical arguments presented above to 2-species systems |155j . One has 
to study two perturbations in each conserved density which due to the inter- 
action are forced to move with the same velocity. Generalizing the analysis 
of the dynamic structure function to two conservation laws one finds that 
the main difference to the case of one conservation law is the evolution of 
two distinct pairs of perturbations out of a single pair. Each pair moves with 
collective velocities given by the eigenvalues of the Jacobian 

/ _d_jA \ 

-D = ^g" ,B • (3.23) 

Corresponding to the two pairs of perturbations a single shock splits into 
two separate shocks, leaving the system in stationary regimes of three distinct 
densities, viz. the left and right shock densities p^2^ enforced by the initial 
state and selforganized intermediate shock densities p"^'^. The equality of 
each pair of shock velocities (given by the general expression ()3.19p applied 
to each single shock) and the requirement that the velocity of the left shock 

must be smaller than the velocity of the right shock determines the 
intermediate densities p^'^. Requiring that all perturbations be absorbed in 
the shock one arrives at the condition 

vt{pi)>v':>v-{p), vt{p)>v'^>vt{p2) (3.24) 

for shock stability in driven diffusive systems with two conservation laws. 
Violating one of these conditions leads to rarefaction waves which have par- 
tially been described [155J . The general features discussed here are confirmed 
by Monte-Carlo simulation of a two-lane model with a conserved density on 



29 



each lane related to the models of Refs. [SH EO], but with periodic bound- 
ary conditions and different choice of hopping rates respectively. A complete 
description of the evolution for all possible initial states has not yet been 
achieved. 

4 Critical Phenomena 

It is well-known that in thermal equilibrium one-dimensional systems with 
finite local state space and short range interactions do not exhibit phase tran- 
sitions at positive temperatures, only at T = long range order may exist. 
From a dynamical viewpoint there are no thermal fluctuations at T = in 
a classical system. In a terms of a stochastic process that means that all 
transition rates are zero. Conversely, if a transition rate is non-zero, some 
dynamics - not necessarily satisfying detailed balance - is going on and it 
has been conjectured that quite generally a system with strictly positive 
transitions rates and local interactions can have at most one stationary dis- 
tribution, which is often rephrased by saying that there can be no phase 
transition in a one-dimensional system with strictly positive rates. One has 
in mind an infinite system since in a finite system dynamics with strictly 
positive rates are always ergodic and the conjecture is trivially true. 

To rationalize the conjecture one imagines, in the simplest case, two po- 
tentially stationary distributions characterized by a different value of the 
order parameter. An example is the Ising model where the order parameter 
is the magnetization, which can take two different values below the critical 
temperature in two or higher dimensions. The reasoning behind the positive 
rates conjecture is the difficulty to imagine a local mechanism that eliminates 
islands of the minority phase (created constantly by thermal fluctuations in 
a region where the other phase dominates) since in one dimension energetic 
effects due to line tension play no role. A local mechanism cannot detect the 
size of a minority island, therefore such an island can grow indefinitely and 
destroy the majority phase. Since noise (implied by strictly positive rates) 
can always create such islands there seems to be no possibility to keep the 
majority phase stable against fluctuations. In a certain "natural" class of 
systems with nearest-neighbour interaction this conjecture has been proved 
rigorously some time ago |156j . 

Therefore it came as a surprise that P. Gacs constructed a model on 
the infinite lattice which violates the positive rates conjecture |157t I158j . 
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However, both the model and the proof that there is a phase transition is 
rather comphcated |159j . requiring either a very large local state space or 
a very large interaction range, and the quest for simple models with this 
property continues to stimulate research. As a guideline we note that the 
conjecture is clearly true for dynamics satisfying detailed balance with respect 
to a local interaction energy as in this case the stationary distribution is just 
the usual equilibrium distribution and the argument underlying the positive 
rates conjecture applies. Hence one should look for models that either violate 
detailed balance or have a nonlocal interaction energy, but local dynamics. 

We address the question of phase transitions in a broader sense by ask- 
ing whether phenomena associated with phase transitions such as divergent 
length scales or spontaneous symmetry breaking may occur far from equilib- 
rium. Divergence of some correlation length does not necessarily require the 
existence of more than one stationary distribution for a given set of system 
parameters. One the other hand, one could have a parameter range with 
two or more stationary distributions which are not related by any symmetry, 
but the transition into this regime would be associated by a divergent length 
scale. 

On the other hand we restrict ourselves to driven diffusive systems with 
one or two conservation laws, thus skirting the issue of phase transitions 
in non-conservative models, addressed in |103t I160t 1161] , and also avoiding 
systems kept out of equilibrium without a having current in the conserved 
density. This could be achieved e.g. by coupling symmetric hopping dy- 
namics to heat baths of different temperatures, a scenario not envisaged in 
the context of driven diffusive systems. Also non-conservative processes with 
absorbing states such as the contact process [TH El Hi] fall in the class of 
systems not considered here. We refer to models of this type only where it 
serves to illuminate the properties of closely related driven systems. 

By definition, conservative systems have a continuum of stationary states 
(characterized by the value of the order parameter) and hence the critical 
phenomena we review concern transitions between different stationary dis- 
tribution at the same value of the order parameter and coexistence of macro- 
scopic stationary domains where the order parameter takes different values. 
The domain walls separating these domains are the shocks discussed in the 
previous section. Hence the stability of domain walls is intimately connected 
with the existence of phase separation. 
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4.1 One conservation law 



In a system with a conserved density the positive rates conjecture does not 
apply by definition, as transitions violating density conservation have zero 
rate. Hence there seems to be no reason to pursue the question of existence of 
phase transitions in conservative systems. Moreover, the discussion of shock 
stability presented above shows that a stable domain wall separating regions 
of different value of the order parameter may exist, thus actually suggesting 
the possibility of phase coexistence and long-range order. However, by the 
same reasoning it is clear that in a system with generic current density re- 
lation one cannot have two domain walls which would be necessary to have 
macroscopic phase separation in a translation invariant system: If say, an 
upward shock from density pi to density p2 at the left boundary of the re- 
gion of higher density p2 is stable by the criterion ()3.22|) . then the downward 
shock p2 to pi at the other boundary of the region of higher density would be 
unstable by virtue of the same criterion. Thus one is forced to conclude that 
phase coexistence in a system with one conservation law may exist only in 
the absence of translation invariance. This is indeed known in systems with 
open boundaries where external particle reservoirs enforce regions of constant 
boundary densities pi,2, separated by a single stable domain wall [201 CHI EH] ■ 
Similarly in a periodic, but not translation-invariant system with a defect 
one may have phase separation since in this set-up the defect may stabilize 
the intrinsically unstable shock. Hence, stability of a shock in a system with 
one conservation law does not constitute a violation of the zero-rate conjec- 
ture for a translation invariant system. Indeed, it has been suggested that 
one-dimensional driven diffusive systems do not exhibit long-range order in 
their steady states 

Yet, several translation invariant models with one conservation law and 
short-range interactions which exhibit a robust phase transition were dis- 
covered. A simple, but non-generic example is a growth model (which can 
be mapped to driven diffusive system according to the strategy explained in 
Sec. 2) where a roughening transition from a phase with a smooth interface 
to a phase with a rough interface occurs |lti2t . This class of models is 
non-generic in so far as there is an intrinsic maximal growth velocity of the 
interface, enforced by a discrete-time parallel update. The smooth phase and 
hence the phase transition disappears if the limit of continuous time is taken 
in these models. 

Addressing the possibility of a roughening transition in systems with con- 
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tinuous time evolution Alon et al |164t I165j proposed a two-species ASEP 
with rates 



Dao = Dob = (1 - g)/2, Dqa = Dbo = Dba = q 



(4.1) 



and annihilation/creation rates 



Doo = Dba = q, Dab = 1- q 



(4.2) 



for the transitions 00 AB and BA — > 00 respectively. These dynamics 
lead to a single conserved "density" S = — N^. In the mapping to a 
growth model A [B) represent a local slope 1 (—1) and represents slope 
0. In order to ensure periodic boundary conditions also in the interface 
representation the model has been studied for S = 0. For small q < qc there 
is a smooth phase where a local mechanism eliminates islands in a fiat region, 
since islands are formed with boundaries that are biased to move towards 
each other. This mechanism applies for islands of all sizes (except completed 
layers) and hence leads to a smooth interface, consisting mainly of vacancies 
in the lattice gas picture. Above Qc ~ 0.189 the creation of new islands 
overcompensates the disappearance due to their intrinsic tendency to shrink 
and the system is in a growing rough phase with a finite fraction of particles. 
The growth model is in the universality class of the KPZ equation which 
is represented by the standard ASEP. The critical behaviour at qc which is 
related to directed percolation is discussed in detail in |lti41 lltiSj where also 
a version of the model without constraint on the local slope (corresponding 
to the absence of exclusion in lattice gas language) is discussed. There is 
spontaneous symmetry breaking in the smooth phase which can be quantified 
by introducing either a colouring scheme, giving vacancy clusters between 
A, B pairs a colour, or by introducing a nonconserved order parameter M that 
makes use of the interface representation. The critical exponent 6 associated 
with the vanishing of the order parameter by approaching the critical point 
from below. 



is a new exponent, not yet understood in the framework of directed percola- 
tion. 

For this transition to exist it is crucial that no chipping of particles from 
the interface in a locally fiat environment may occur. This is ensured by 
imposing a kinetic constraint on the local dynamics, viz. setting the rate 




(4.3) 
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Dqq for the process 00 — > BA to zero. However, this process does not violate 
the single conservation law. Hence this model does not provide a coun- 
terexample against the zero-rate conjecture, as applied to systems with one 
conservation law. By regarding the various local transitions of this strongly 
nonequilibrium process as induced by thermal activation from heat bathes at 
different temperatures one is led to conclude that the local mechanism that 
guarantees bounded growth of regions of the minority phase is bought with 
a zero-temperature condition on the chipping process. 

In the presence of chipping with a rate Dqo = p the interface attains a 
negative stationary growth velocity for some value of q that depends on p. 
An interesting phenomenon then occurs if the dynamics of the interface is 
constrained by the minimal height condition that hi > for all times and 
all lattice points i |166j . For negative growth velocity the interface is driven 
towards the hard wall located at the height level hi = which the interface 
cannot penetrate. In the special case of p = 1 — q > 1/2 the model satisfies 
detailed balance with respect to the energy 

E = j:h, (4.4) 

i=l 

which is the area under the interface. The stationary probability of finding 
an interface configuration h = {hi, . . . ,hi) is given by 

P{h) = (g/(l - q)f^^^/ZL (4.5) 

with the partition function Zl = J2config{l/ At g = 1/2 the interface 
has mean velocity zero, for g > 1/2 the interface grows. Hence the expression 
(j4.5j) diverges in time and becomes meaningless as stationary distribution 
for g > 1/2. In the growth regime the interface roughens, with a dynamical 
behaviour in the ubiquitous KPZ universality class. For q < 1/2 the interface 
is bound to the hard wall and hence is smooth. Close to the transition point 
exact analysis of the partition function yields an occupation density a of the 
bottom layer h = and a width w diverging as |16(ij 

(^-^{qc-q)\ w~(g,-g)i/3 (4.6) 

The unbinding of the interface at gc = 1/2 is analogous to a wetting transi- 
tion. We remark that this model has a zero-rate constraint by not allowing 
the interface to penetrate the bottom layer. 
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The physics described by these two models can be captured in a general- 
ized KLS-model (|2.11|) . In this model spontaneous symmetry breaking due 
to absence of chipping, and conservation of minimal height with the resulting 
wetting transition can be studied without reference to the height variable. 

Koduvely and Dhar considered the symmetric KLS-model with rates 
D^^ ■= D^^ = Df^ for the hopping event XAOY ^ XOAY ^U7\. The 
analog of the chipping rate in the two-species model is the hopping rate D^*^. 
Setting D^'^ = automatically leads to conservation of minimal height, or, 
more precisely, conservation of the height level of a completed layer. It is 
not necessary to stop the dynamics by a separate rule involving the local 
height. In contrast to the previous model of Hinrichsen et al, however, the 
interface always remains anchored to the minimal height level at some ran- 
dom position. Hence the properties of the wetting (unbinding) transition 
are described by this model only below the critical point in the dry (bound) 
state. 

Careful numerical analysis of the symmetric model |167j indicates subdif- 
fusive critical dynamics of the unbiased interface with a dynamical exponent 
z ~ 2.5, as opposed to the usual dynamical exponent z = 2 of the unbiased 
Edwards- Wilkinson interface modelled by the symmetric KLS model for any 
D^^ > 0. The physics of the biased interface where studied by Helbing et al 
|168j by considering asymmetric rates 

Dr = r, D^^ = q, Dr=P, Dr = 0. (4.7) 

with p = 1 — g — r. For r = 0, the system satisfies detailed balance with 
respect to the measure (j4.5p . Below the wetting transition q < 1/2 one has 
the exact exponents (j4.5|) . At the critical point the interface is rough and 
one expects a dynamical critical exponent as measured in the related model 
of Koduvely and Dhar. Above the critical point the interface would grow, 
but cannot detach from the minimal height level due to anchoring. Hence 
the measure (j4.5|) is stationary for all q. The stationary interface has a cusp, 
the anchoring point is random. This is an example of spontaneous breaking 
of translational invariance. In the particle picture the steady state is a shock 
measure with extremal limiting densities pi = e, p2 = 1 — e with a sharp 
downward shock at some random lattice site k. In a finite system the small 
quantity e and hence the particle current is exponentially small in system 
size L. The random anchoring point k moves with a speed also exponentially 
small in system size. This structure describes a strongly phase-separated 
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state with an essentially empty region and an essentially full region. For this 
to be valid it is not necessary to require a total average density p = 1/2. It is 
easy to understand this steady state directly from the microscopic dynamics. 
The right edge of the occupied domain, i.e., the right-hand shock is stable 
because a hopping of the rightmost particle is exceedingly unlikely since the 
transition AAO —>■ AO A is forbidden and the configuration OAO where hopping 
is allowed is exceedingly unlikely for the rightmost particle (exponentially 
small in system size). On the other hand, at the left boundary of the domain 
the system essentially behaves like an ordinary ASEP which has a stable 
upward shock in the direction of motion. 

For r > the interface is not anchored anymore, but the minimal height 
condition for any completed layer is still conserved. Along the line r = q one 
observes a transition from a rough growing interface (finite particle current 
J > for q > Qc ^ 0.1515) with KPZ dynamics to a smooth interface with 
spontaneously broken symmetry. The symmetry breaking can be quantified 
in terms of the non-conserved order parameter 



This quantity measures the difference of sublattice densities between the even 
and odd sublattice respectively. In the language of spin systems this is the 
staggered magnetization, playing the role of the order parameter for antifer- 
romagnetic systems. Below the critical point particles either accumulate on 
the even or odd sublattice respectively. Both happens with equal probability, 
but a transition between both kinds of configurations occurs on a time scale 
that diverges exponentially in system size. This is the signature of sponta- 
neous symmetry breaking in a finite system. An explicit calculation of the 
exponentially large transition time is possible in the vicinity of the line g = 
|169j . No reference to the height variable is necessary for measuring M. Nu- 
merical investigation of the model shows that it is in the same universality 
class as the growth model of Alon et al. |168j . 

In neither of the models discussed above the positive-rates condition is 
satisfied. Only kinetic constraints imposed by vanishing rates, analogous to 
the zero-temperature condition on phase transitions in equilibrium, may lead 
to a vanishing current in which case the stability argument for domain walls 
does not apply and phase transitions can occur. Below we shall present an in- 
dependent argument that suggests conditions under which kinetic constraints 
lead to phase separation. Hence so far there is no known simple model with 




(4.8) 



k 
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a single conservation law that violates the positive rates conjecture. 

4.2 Phase separation in two-species ASEP's 

The exact and numerical analysis of steady states of one-species systems 
reviewed above has revealed that phase separation in systems defined on a 
finite ring or on Z may occur if one or more of the following conditions are 
satisfied: 

• (I) there are spatially localized defects reducing the mobility of particles 

• (II) single particles of a different species act as mobile blockages 

• (III) the dynamics have kinetic constraints arising from a nonequilib- 
rium zero-temperature condition. 

The last condition leads to strong phase separation in the sense that one 
domain is fully occupied whereas the other domain is entirely empty. The 
current in the phase separated state vanishes exponentially in the size of 
the particle domain, the separated state exists at any total particle den- 
sity. Conditions (I) and (II) may lead to strong phase separation, but allow 
also for a soft phase separation between domains of different densities. This 
phenomenon sets in only for densities above some critical density pc- The 
steady-state current is nonvanishing and independent of p in the phase sep- 
arated state: Increasing the density leads to an increase of the size of the 
high-density domain, but not to a change of the current. In analogy to 
Bose-Einstein condensation we call the high density domain a condensate, 
the transition at pc is referred to as condensation transition. Notice that 
this characterization refers to the thermodynamic limit L oo. In a fi- 
nite system there is either a current exponentially small in system size (case 
A, strong phase separation) or one has finite-size corrections to the finite 
bulk current (case B, soft phase separation). We remind the reader that soft 
phase separation may disappear above a critical density pc and also for a 
finite density of blocking particles 95 . Strong phase separation is accompa- 
nied by spontaneous breaking of translational invariance, except if caused by 
condition (I) where translational invariance is explicitly broken. 

Strong phase separation has been found also in homogeneous systems on 
a ring where neither of the conditions (I) - (III) is satisfied, but where there 
is a second species of particles with finite density (HOI I143t ll7Uj . Hence we 
add a further sufficient condition for the possibility of phase separation 
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• (IV) the system has two or more conservation laws 

We remark that all the conditions (I) - (IV) in some way or other impose 
local constraints on the dynamics of the driven diffusive system. This appears 
to be a general requirement for phase separation in generic driven diffusive 
systems. The size of the local state space and the range of interaction appear 
to be irrelevant if one of the conditions (I) - (IV) is satisfied. 

Using a four-states model which is equivalent to a two-lane model with two 
conserved densities Lahiri and Ramaswamy [301 1171] address the question of 
phase separation in terms of the stability of crystals moving steadily through 
a dissipative medium, e.g., a sedimenting colloidal crystal. In a certain limit 
(large particle radius or small elastic modulus of the suspension) experiments 
suggest instability of such a crystal. Numerical analysis of the lattice model, 
however, reveals a transition to a stable regime, corresponding to strong 
phase separation. In the simpler 2-species ASEP (j3.9|) or in a more symmetric 
model with rates [1701 |I72j 

Dao = Dqb = Dba = 1, = Dbo = Dab = q (4.9) 

the mechanism for strong phase separation for g < 1 is very transparent. Here 
strong phase separation refers into separation of three pure macroscopic do- 
mains, each consisting of essentially only one particle species or empty sites. 
For simplicity we assume = , but this is not necessary for the phe- 
nomenon to occur. Prepare a phase-separated block which we symbolically 
represent by . . . OOOAAAAAABBBBBBOOO .... One observes the following: 
(i) The 0\A interface is stable by the criterion ()3.22|) since due to the absence 
of i?-particles one has the dynamics of the usual ASEP (with a bias to the 
right) in the vicinity of this domain wall, (ii) The B\0 interface is stable 
for exactly the same reason {B particles have a bias to the left) (iii) The 
A\B interface is stable since in the absence of vacancies S-particles act like 
vacancies w.r.t. the local dynamics of the A-particles and vice versa, (iv) 
Since each domain wall is stable (only small fluctuations extended over a 
finite range of lattice sites evolve at the phase boundaries) the assumption 
used in the argument remains valid for all times. 

It is clear that this model can be extended to an arbitrary number of 
conserved species and does not require equal density for each particle species. 
However, for A^"^ = A^-^ = L/3 the dynamics can be shown to satisfy detailed 
balance w.r.t. an equilibrium measure with a non-local interaction energy of 
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The corresponding partition is proportional to system size L, rather than 
exp(/L), since fluctuations occur only in a finitely extended region around 
the three domain walls. The "temperature" associated with equilibrium mea- 
sure is given by kT = — l/(lng). It diverges at g = 1 which corresponds to 
the disordered state of symmetric diffusion of first- and second-class parti- 
cles. The analysis of the model for g > 1 is similar, with the role of A and 
-B-particles interchanged. The behaviour of the model close to the transition 
point q = 1 has been investigated |173j and some critical exponents have 
been determined numerically. 

The related but distinct model ()3.9|) of Arndt et al. has strong phase 
separation for g < 1, with essentially nonfiuctuating 0\A and B\0 domain 
walls, and an y4|i?-interface similar to the one in model ()4.9|) . For g > 1 the 
behaviour of the model is more intricate. Numerical and mean-field analysis 
|174j suggests the existence of soft phase separation up to a critical value 
Qc = I + 4:Xp/ (1 + 2p). There is a condensate of density 1, but consisting of 
both species of particles, but vacancies. The other "fiuid" phase has density 
< 1, with particles of both species and vacancies distributed apparently 
similarly to the disordered phase for g > gc. Notice that inside the condensate 
particles also fiow, but with the dynamics of the usual ASEP, as the B 
particles act like vacancies in the usual ASEP. The number of A and B 
particles is on average equal and the condensate essentially behaves like the 
usual ASEP with open boundaries in the maximal current phase, except 
that the system size M of this ASEP corresponds to the slightly fiuctuating 
cluster size. The stationary current for both particle species is non-zero and 
approximately given by the value J = (g — l)/4 expected from the ASEP in 
the maximal current phase. The position of the condensate fiuctuates on the 
lattice. For an unequal average density of A and 5-particles, the scenario 
as described here remains essentially unchanged, except that the condensate 
has a finite drift velocity ,175j . 

The exact stationary distribution of the model (|3.9p can be calculated 
using the matrix product ansatz. For the grandcanonical ensemble with equal 
densities of both particle species it has been shown |176j that the apparent 
condensation transition is a crossover effect. For sufficiently large lattice 
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one would observe a distribution of clusters, but not a single macroscopic 
condensate. Using the parameter a = (1 — q^^ — A)/A and the fugacity ^ 
controlling the density the exact current J and the density p are given by 

"^^^^ " 1 + a2 + 2a(l + a2)e - (1 - a^)^! + 



P(0 = . r. -— — — ^=^.(4.12) 



a(l + a)C, 


(l + a)v/l + 4a^- (1 - a) 




VI + 4a^ 


1 + a2 


+ 2a(l + a2)^ - (1 - a^) ^1 + 4a^ 



This expression has a very remarkable property: Inside the apparently con- 
densed phase at g = —a = 10/9 the derivative J'{p) which enters the Jacobian 
(|3.23p and hence determines the collective velocities of the two-species sys- 
tem has a change of order 1 arising from a change of order 10~^^. One would 
need a lattice of the order of 10™ sites to actually observe the breakdown of 
the condensation and see the full distribution of clusters of various lengths. 

Clearly, a crossover scale of this magnitude is of no relevance for the oc- 
currence of soft phase separation in a real system. Any finite sample would 
exhibit a phase separated state. However, the huge number 10™ character- 
izing the typical scale of the cluster size distribution is specific for the model 
()3.9p . As shown in |177j the crossover scale depends sensitively on nonuniver- 
sal parameters which are tunable in some models |178j . Hence other models 
may have parameter ranges with a mean cluster sizes of the order, say, 10^. 
Such a crossover scale would render numerical results from computer simula- 
tions of the steady state for a real quasi one-dimensional system ambiguous: 
A computer simulation of a realistic model for a real system with less than 
10^ particles could predict soft phase separation, whereas actual experiments 
done on a macroscopic sample with more than 10^ particles could yield the 
contradictory result that there is only a disordered phase. Yet, one could not 
conclude from this observation that the model is inappropriate to describe 
the real system since in a smaller experimental sample of the same system 
computer simulations and experimental observations may agree. 

We remark that a technical assumption in the exact calculation of Rajew- 
sky et al |175j has been proved in [TTS] . The validity of the result also in the 
canonical ensemble has been challenged |175j . even though exact analysis of 
the fluctuations in the particle density strongly suggest that there is no true 
condensation in the canonical ensemble |176j . In another model introduced 
by Korniss et al a two lane extension of a 3-species driven system was stud- 
ied jl801 1181] . It has been suggested that while for this model the one lane 
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system does not exhibit phase separation |182j . this phenomenon does exist 
in the two lane model. The studies rely on numerical simulations of systems 
of length up to 10"^. However, no theoretical insight is available as to why 
phase separation in this two-lane model should persist in the thermodynamic 
limit. 

To conclude these general considerations we note that in many real quasi 
one-dimensional systems the particle number is in the range 10 . . . 10^. Sim- 
plified models of such systems are accessible to numerical simulations of the 
steady state. This raises the further question to which extent phase sep- 
aration in finite systems is obscured by a too large intrinsic width of the 
domain wall separating the condensed domain from the "fluid" low-density 
domain. While quantitative predictions for the relaxation modes resulting 
from the coarse-grained domain wall theory in single-species systems with 
open boundaries have been verified numerically for the ASEP on small lat- 
tices of only C(10) sites jl83j there is no systematic finite-size study of soft 
phase separation. 

4.3 Criterion for phase separation 

In view of the exact analysis of the model ()3.9p it is clear that numerical evi- 
dence for soft phase separation may be rather subtle and indeed be mislead- 
ing. It would thus be of great importance to find other criteria, which could 
distinguish between models supporting phase separation from those which 
do not. Phase separation is usually accompanied by a coarsening process in 
which small domains of, say, the high density phase coalesce, eventually lead- 
ing to macroscopic phase separation. This process takes place as domains 
exchange particles through their currents. When smaller domains exchange 
particles with the environment with faster rates than larger domains, a coars- 
ening process is expected, which may lead to phase separation. 

An approach that quantifies this mechanism and yields a criterion for 
phase separation in terms of the current leaving the domains is proposed by 
Kafri et al jl84j . The current-criterion is readily applicable even in cases 
which cannot be decided by direct numerical simulations. In order to explic- 
itly state the criterion one distinguishes systems with a vanishing current of 
a finite domain of size n 

Jn^O (case A) (4.13) 
from systems finite-size corrections to a finite asymptotic domain current Joo 
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of the form 

Jn = ^oo(l + h/n") (case B). (4.14) 

to leading order in 1/n. For simplicity we assume here domains with van- 
ishing drift velocity in which case the current inside the domains equals the 
outgoing current. More generally one has to distinguish the two currents 
leaving the cluster at the right and left boundary respectively. 

For 6 > the current of long domains is smaller than that of short ones, 
which leads to a tendency of the longer domains to grow at the expense of 
smaller ones. The current criterion asserts that phase separation exists only 
in the following cases 



Jn for n ^ oo (case A) (4-15) 
Jn Joo > (case B) (4.16) 

for either a < 1 and 6 > or a = 1 and b > 2. In case A one has strong phase 
separation for any density, whereas in case B one has soft phase separation 
at any density for a < 1 and above a critical density 

P==^ (4.17) 

for cr = 1. The fluid regime has particles with density pc- Hence in a 
finite system the macroscopic size of the condensate in the phase-separated 
regime is determined by the system parameter b. For an asymptotic decay 
faster than 2/n there is no condensed phase, the system is disordered for all 
densities. 

Models with two conservation laws for which J„ decays exponentially to 
zero with n (case A) have been reviewed above and indeed were shown to 
exhibit strong phase separation at any density. In the model (j4.7|) with a 
single conservation law the current of particles out of the left domain wall, 
i.e., the current opposite to the bias of the individual particles, decays expo- 
nentially with domain size for r = 0, as demonstrated for the usual ASEP 
with appropriately chosen open boundary conditions ^E^- The domain size 
dependence of the current flowing away from the right edge of a cluster (in 
direction of the bias) is to leading order not a self-organized quantity, it is de- 
termined strongly by the interaction of the particle at the edge of the cluster 
with the surrounding particles. The interaction range is one lattice site and 
hence the current J„ out of the right edge of a cluster of density 1 becomes 
to leading order independent of the cluster size for n > 2. It vanishes due to 
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the kinetic constraint r = (zero-temperature condition). According to the 
criterion one expects strong phase separation, in agreement with the result 
reviewed above. 

For Joo 7^ (case B) we note that in a system with two conservation 
laws the current inside a cluster organizes itself to a value determined by 
the dynamics of the reduced system with only one conservation law resulting 
from the absence of vacancies. This reduced system has open boundaries with 
in- and outflow of particles such that the system is in the generic maximal 
current phase of the reduced system. It is assumed that the current flowing 
through a block is given by its steady-state value and is independent of its 
neighboring blocks. This may be justified by the fact that the coarsening time 
of large domains is very long, and the domains have a chance to equilibrate 
long before they coarsen. 

In case B one expects generically a = 1 for the following reason: (a) In 
a periodic system the leading finite-size corrections to the current Joo in a 
canonical ensemble is given by J„ — Joo = ^•J!L'^/i'^^) |186t 1178] . Here J^ 
is the curvature of the current-density relation and k = {{N"^) — {N)'^)/L 
is the nonequilibrium analog of the thermodynamic compressibility which is 
assumed to be finite, i.e., one assumes sufficiently rapidly decaying correla- 
tions as was implied above in the derivation of the collective velocity which 
also requires finite compressibility, (b) There is a universal ratio c* of the 
finite-size corrections to the current in the maximal current phase of a driven 
diffusive system (which describe the dynamics inside the growing domains) 
and the finite-size corrections of the canonical ensemble of a periodic system 
[r§7j . This yields leading finite-size corrections of the form ()4.13p with a 
parameter b entirely determined by the universal constant c* and the macro- 
scopic quantities J^ and k. The value of c* = 3/2 has been obtained from 
the exact solution of the ASEP with open boundaries [T^ l2Uj . 

We stress that by definition 6 is a quantity that itself does not depend on 
system size. For systems with unknown stationary distribution the reduced 
dynamics inside a cluster allows for a simple numerical measurement of b by 
studying the finite-size corrections of the stationary current in the reduced 
open system of length n. One neither needs huge lattices nor is one faced 
with the problem of slow relaxation of the phase separation process in the 
full system. Applying the criterion to the model ()3.14|) yields the exact value 
6 = 3/2 and hence one expects no condensation, in agreement with the exact 
result. For the two-lane model of Korniss et al |180j one obtains numeri- 
cally b ^ 0.8 |184j and therefore one expects no condensation in contrast to 
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the results of the Monte-Carlo simulation of the full model with 10'* lattice 
sites. A three-states model with KLS-dynamics (|2.11|) inside the clusters has 
been shown to have b > 2 |178j which suggests the existence of soft phase 
separation in driven diffusive systems with two conservation laws. 

The criterion presented above emerges from a careful analysis of the zero- 
range process (ZRP) which could be viewed as a generic model for domain 
dynamics in one-dimension |184j . 

Depending on the rates Wn the model may or may not exhibit conden- 
sation in the thermodynamic limit, whereby the occupation number of one 
of the boxes becomes macroscopically large. Clearly the rate Wn must be 
a decreasing function of n in order for larger blocks to be favored and to 
support condensation. It is known jHl llU6j that condensation occurs at any 
density when w„ — >■ with n — cxo, or when it decreases to a non- vanishing 
asymptotic value as b/n" with a < 1; no phase separation takes place for 
0" > 1; for (7 = 1 phase separation takes place at high densities only for 
b > 2. This model may be used to gain physical insight into the dynamics 
of driven one-dimensional systems. Occupied boxes represent domains of the 
high density phase. The currents leaving domains are represented by the 
rates of the ZRP. This is done by identifying the rate Wn associated with a 
box containing n balls with the currents Jn leaving a domain of n particles. 
A bias in the currents to a certain direction may be incorporated through 
a bias in the ZRP dynamics. The existence of a box with a macroscopic 
occupation in the ZRP corresponds to phase separation in the driven model. 
The distribution of occupation numbers obtained from the ZRP was shown 
to agree with the domain size distribution of the model ()3.14|) |184j . 

It is remarkable that extending the asymptotic expansion ()4.14|) by a 
(nonuniversal) next-leading term c/n"^ leads to an extremely sensitive depen- 
dence of the mean domain size ^ on c \177\ . The quantity ^ exhibits a sharp 
increase of a few orders of magnitude over a narrow range of values of c. This 
reflects itself in large (but finite) blocks and an apparent phase separation in 
direct simulations. 

4.4 Coarsening 

The previous discussion has focussed on the stationary properties of phase 
separation. The dynamics below the critical density are expected to be de- 
scribed on the hydrodynamic scale by the mechanisms reviewed in Sec. 3. 
Above the critical density a natural set-up is to start with particles uniformly 
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distributed at the supercritical density p > pc- In the beginning the excess 
particles condense at a few random sites. Thus there are several clusters 
which are essentially immobile. On the remaining sites, the distribution re- 
laxes to its critical stationary distribution with p = pc- With increasing 
time the larger clusters will gain particles at the expense of the smaller ones, 
causing the clusters to merge. Eventually only a single cluster containing 
all excess particles survives, which is typical for the stationary distribution 
(Fig. 6). 

In case A (strong phase separation) the evolution of this state will pro- 
ceed by slow diffusion against the bias, in which, for example an A-particle 
traverses the adjacent domain of 5-particles. Using standard mean first 
passage time calculations the time t necessary for penetrating the complete 
5-domain of size n can be shown to be of order where g > 1 is the hopping 
asymmetry, i.e., the ratio of hopping rates inside the i?-domain. After that 
the A-particle travels with finite average velocity to the edge of the next A 
domain. This happens after a time of order L which can be neglected com- 
pared to the penetration time. Therefore the coarsening time is of the order 

and inverting this relation yields a logarithmic growth law for the mean 
domain size in the strong phase separation \112\ 

n(t) ~ Int/lng. (4.18) 

This growth law is valid also in higher dimensions |188j . 

The zero-range picture may be used in order to study the coarsening 
dynamics of domains in the case of soft phase separation. We set Joo = 1 
which only fixes an uninteresting microscopic time unit for coarsening. In 
[189 the coarsening dynamics are described by studying with random walk 
arguments the loss and gain of particles in neighbouring clusters, mediated by 
stationary transport in the fluid phase in between. The time that particles 
lost from a cluster spend in between clusters is of the same order as the 
time required to dissolve a cluster completely. Hence this time scale is the 
appropriate coarsening time scale. In the totally asymmetric case excess 
particles leave a cluster site with n particles at a rate h/n and therefore the 
typical required to dissolve such a cluster is given by tain) ~ /h. Inverting 
this scaling relation yields the growth law 

n{t) ~ {htf^ (4.19) 

for the mean cluster size n. From ()4.19|) one reads off the dynamical exponent 
z = 2 m the asymmetric condensed phase. 
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Because of the recurrence of a 1-d random walk in the case of symmetric 
hopping, a particle that got lost from a cluster is very likely to return to the 
cluster it left, in contrast to the asymmetric case. It will be trapped by the 
next cluster only with a probability inversely proportional to the diffusion 
distance, i.e., of the order (p — p^jra. So the typical time of a particle to 
leave a cluster is of the order irP'/{h{p — Pc))- Hence the coarsening time 
scale is given by tgin) ~ ta{p — Pc)/m ~ m^/[b{p — pc)]. This results in the 
growth law 

n(t)~[6(p-p,)t]i/3 (4.20) 

and dynamical exponent z = 3 in the symmetric case. A similar growth law 
which can be rationalized using analogous arguments is known for equilibrium 
Kawasaki dynamics (j2.12|) at very low temperatures |19U1 1191j . For biased 
coarsening dynamics with conserved order parameter one has z = 2 |192H193] . 
in agreement with ()4.19|) . The coarsening stops when only one macroscopic 
cluster is left. The typical time of a macroscopic fluctuation of the cluster 
size diverges exponentially with the system size L |189j . 

The growth laws (|4.19j) . (j4.2Uj) were confirmed by Monte-Carlo simula- 
tions of the ZRP |189j . The critical exponents were obtained independently 
by numerically studying the second moment of the local density |194j . Also 
a universal scaling function for the cluster size distribution was obtained. At 
the critical density itself Monte Carlo simulations suggest dynamical expo- 
nents Zc ^ 3 for the asymmetric case and ~ 5 for the symmetric case. A 
theoretical derivation of these exponents is lacking. 

5 Conclusions and open questions 

During the last decade the study of one-dimensional driven diffusive systems 
has contributed significantly to the understanding of critical phenomena far 
from equilibrium, where "understanding" not only refers to the characteriza- 
tion of nonequilibrium universality classes in terms of critical exponents, but 
also to the identification of some of the coarse-grained dynamical mechanisms 
that generate these critical phenomena. Exact results for simple model sys- 
tems such as the ASEP have played a crucial role in advancing and shaping 
this understanding. They have provided deep and detailed insights in robust 
generic phenomena which could then be generalized to more complicated 
systems. 
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This includes the derivation of exact hydrodynamic equations and asymp- 
totic coarsening laws from the microscopic laws of interaction, both rigorously 
and using more intuitive physical arguments based on the application of con- 
servation laws and using random walk arguments. Thus complicated col- 
lective phenomena such as shocks can be viewed as effective single-particle 
excitations with simple properties. Universal fluctuations determining the 
dynamics on intermediate scales between the microscopic lattice scale and 
the macroscopic Euler scale have been probed using test particles: tracer 
particles, second class particles, impurity particles. Thus some similarities 
to the still unresolved problem of localized blockages could be established. 
The occurrence of stationary critical phenomena, viz. spontaneous symme- 
try breaking, long-range order and phase separation of two types, soft and 
strong could be linked to microscopic properties of the dynamics, listed in 
conditions I - IV of Sec. 4. Stability criteria ()3.22|) . ()3.24|) for microscopi- 
cally sharp domain walls and coarsening dynamics of domains ()4.15|) . ()4.16p 
as well as the phase diagram of open one-species systems ()1.5|) can be phrased 
directly in terms of the macroscopic current. We conclude that the current- 
density relation (which is determined by microscopic parameters) is a central 
quantity determining the large scale physics of driven diffusive systems. This 
is fortunate in the sense that the current is usually relatively easy to measure 
or calculate. It is crucial to know the current exactly, mean-field approxima- 
tions are not likely to produce even qualitatively (let alone quantitatively) 
correct features of the system, unless by chance certain analytic properties 
of the mean field current happen to coincide with those of the exact current. 
However, this is not guaranteed even if bulk correlations in the stationary 
state are short ranged. 

A major open question concerns the precise relationship of the various 
current criteria and the microscopic criteria for phase separation. Addressing 
this issue leads to series of interrelated problems, each of which is interesting 
in its own right. 

It is not clear how to deduce the stability of a domain wall in the sta- 
tionary condensed state or during the coarsening process (where the current 
is also stationary and independent of the density) from the stability crite- 
ria ()3.22|) . ()3.24|) . This requires a more careful analysis of the behaviour of 
fluctuations in systems with two conservation laws and is intimately linked 
to a hydrodynamic description of the dynamics above the critical density, 
which is an open problem even within the framework of an effective ZRP 
description. The hydrodynamic treatment of the blockage problem in the 
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usual ASEP |149j may provide some insight. 

Moreover it would be highly desirable to have a unified picture which 
allows for an application of the current criteria ()4.15p . ()4.16p for phase sep- 
aration to single-species systems satisfying one of the conditions (I) - (III) 
and thus to predict from the dependence of the domain current on the mi- 
croscopic parameters at which value of these parameters phase separation 
sets in. This requires a proper definition of the current out of a domain in 
system with one conservation law which might then answer also the question 
whether soft phase separation is possible if neither of the conditions (I) - 
(III) are satisfied. In a single-species system this condensation phenomenon 
is reminiscent of spontaneous traffic jams in automobile traffic fiow. Indeed, 
traffic models with nonconserved internal degrees of freedom are known to 
exhibit soft phase separation |1951 1196| 11971 1198j , but the minimal require- 
ment for the existence of the phenomenon is not entirely clear. The role of 
non-conserved internal degrees of freedom in critical phenomena needs fur- 
ther clarification also in the theory of boundary-induced phase transitions 
which requires some extension |199j . 

In this context it would be interesting to try to predict phase separation 
from the properties of the current, using ()3.22p . ()3.24p . This has not been 
attempted yet even in systems where both the current and the existence of 
phase separation are known and might shed light on the possibility of phase 
separation between domains of nonextremal densities pi ^ 0,1 in translation 
invariant systems. The stability criterion (j3.22j) rules this out for generic 
current-density relations in systems with one conservation law, but it is not 
tantamount to a no-go theorem as a current which is constant in some density 
range may allow for such phase separation. 

A further promising and closely related direction of research concerns the 
hydrodynamics of systems with more than one conservation law. The lack 
of a full hydrodynamic description of particle systems with two conservation 
law still constitutes a major gap in understanding both dynamic and critical 
stationary critical phenomena. The selection of the physical solution using 
a regularization by adding a phenomenological viscosity term is not fully 
understood, as the nature of such a term might not be as arbitrary as for 
systems with one conservation law. This may be of importance not only for 
the bulk critical phenomena reviewed here, but also for boundary-induced 
spontaneous symmetry breaking |2UUj and steady state selection in open 
systems. There is no theory of boundary-induced phase transitions with an 
extremal principle analogous to (jl.5j) that could explain quantitative features 
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of the phase diagram in terms of effective boundary densities as is possible 
for single-species models. Hence it is difficult to make predictions of the 
stationary behaviour in open systems. For systems with more than two 
conservation laws the problem is likely to be even more intricate, but possibly 
also even more fascinating. 

Universal properties of fluctuations which manifest themselves on scales 
below the Euler scale are not readily accessible with the analytical methods 
reviewed in this article. However, with the tools of Bethe ansatz and random 
matrix theory they have become amenable to exact treatment in the frame- 
work of the ASEP POIJ. It is natural to try to apply these techniques 
to the general 2-species ASEP, or at least to the integrable cases. Univer- 
sal quantities - critical exponents and scaling functions, but also universal 
dynamical mechanisms such as evolution of shocks and coarsening - derived 
from studying driven diffusive systems are ultimately expected to be useful 
also in the investigation of real systems such as those listed above. Therefore 
one needs to understand the role of the lattice in the phenomena discussed 
above. Passing to driven diffusive system defined in spatial continuum by 
taking appropriate limits may give an answer. 
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Figure Captions: 



Fig. 1: Asymmetric simple exclusion process and related zero-range 
process. In the ASEP particles on a lattice hop with rates Dr,i to the 
right and left respectively, provided the target site is empty. In the case 
of open boundaries they are created or annihilated with rates a,l3,'~f,6 at 
the boundary sites 1,L as indicated in the figure. Reflecting boundaries 
correspond to a — P — ^ — S — 0, in the case of periodic bound- 
aries particles may hop between sites L,l. In the associated zero-range 
process consecutive particles correspond to sites, the length of empty 
intervals between them to occupation numbers. Hopping of particle k from 
a given in the ASEP corresponds to hopping between sites k — 1, k in the ZRP. 

Fig. 2: Mapping of the ASEP to a restricted solid-on-solid (RSOS) growth 
model. An empty (occupied) site corresponds to a slope +1 (-1) in the 
associated interface. Hopping to the right (left) thus is equivalent to a 
deposition (evaporation) of an atom. The broken line at the bottom of 
the interface marks a completed layer. In the absence of chipping the 
interface cannot shrink below a completed layer which yields a minimal 
height model. For an anchored interface the minimal height of the interface 
equals the height of the top-most completed layer. In two-species mod- 
els the definitions are analogous, but with allowed local interface slopes 0, ±1. 

Fig. 3: Main features density profiles in a periodic system in the presence of 
a localized defect or a mobile impurity at site L. Below the critical density 
there is an exponentially decaying profile behind the disturbance, above 
the critical density one has a macroscopic "traffic jam" extending over a 
finite fraction of the total length L. Depending on the specific system, the 
density profile may have some extra structure (not shown here) around 
the disturbance. The position of the shock at the end of the high-density 
domain fluctuates. 

Fig. 4: Time-evolution of the ASEP on the Euler scale, starting from 
time t — (a) with an unstable shock and a region of positive slope. Any 
region of positive slope evolves into a shock after some time (b) because of 
the flow of localized perturbations (microscopic picture) or characteristics 
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respectively (macroscopic description). The evolution of the unstable shock 
is not uniquely defined by the viscosity-free hydrodynamics. Both evolutions 
(b) and (c) are solutions of the inviscid Burgers equation for initial state (a). 
The physical solution selected by the spatially homogeneous ASEP is shown 
in (b). The time evolution shown in (c) corresponds to the ASEP with a 
defect which is governed by the same homogeneous Burgers equation |149j . 

Fig. 5: Schematic phase diagram with second order phase transition line 
(broken curve) between the disordered phase (A) (growing KPZ-interface) 
and ordered phase (smooth interface) (B) with spontaneously broken 
Z2-symmetry and non- vanishing order parameter M (14. 8|) . At r = 0, there 
is a transition at q = 1/2 with strong phase separation and spontaneous 
breaking of translational invariance (bold line C). 

Fig. 6: Coarsening of domains during phase separation. The regions be- 
tween domains have relaxed to their stationary critical density. Each domain 
of size Ur has an outflowing current jr determined by the dynamics inside 
the domain. In order to have domain growth the current should decrease 
with domain size, as indicated in the figure for domains 1,2,3. Domain 2 is 
expected to be eaten up by domains 1,3 after some time. 
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Figures: 
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